function harmamp = harmamp(fft, pitch, nyquist)

len_fft = length(fft(:,1));

for i=1:length(pitch)
    pitch_frame = round(pitch(i)*len_fft/nyquist);
    for j=1:5
        harmamp(j,i) = fft(j*pitch_frame, i); 
    end
end